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In this paper we calculate the non-perturbative Euler-Heisenberg Lagrangian for massless QED 
in a strong magnetic field H, where the breaking of the chiral symmetry is dynamically catalyzed 
by the external magnetic field via the formation of an electro-positron condensate. This chiral 
condensate leads to the generation of dynamical parameters that have to be found as solutions of 
non-perturbative Schwinger-Dyson equations. Since the electron-positron pairing mechanism leading 
to the breaking of the chiral symmetry is mainly dominated by the contributions from the infrared 
region of momenta much smaller than %/ eH, the magnetic field introduces a dynamical ultraviolet 
CNj ' cutoff in the theory that also enters in the non-perturbative Euler-Heisenberg action. Using this 

action, we show that the system exhibits a significant paraelectricity in the direction parallel to the 
■ magnetic field. The nonperturbative nature of this effect is refiected in the non-analytic dependence 

' of the obtained electric susceptibility on the fine-structure constant. The strong paraelectricity in 

the field direction is linked to the orientation of the electric dipole moments of the pairs that form 
the chiral condensate. The large electric susceptibility can be used to detect the realization of the 
magnetic catalysis of chiral symmetry breaking in physical systems. 

PACS numbers: ll.30.Rd, 12.38.Lg, 81.05.ue 

I. INTRODUCTION 

O ' The study of rclativistic quantum theories in the presence of strong magnetic fields has been an active topic of 
j-^^l research in the physics literature over many years The observation of magnetars' surface magnetic fields in the 
range of lO^^G — lO^^G 0, together with estimates of inner magnetic fields between lO^^G and 10^°G, for star cores 
' made of nuclear Q and quark matter Q respectively, have contributed to renovate the interest in this area. In 
I> . addition, the production of very strong magnetic fields ^ 277i^(~ lO^^G) in non-central heavy- ion collisions at RHIC 
' 0, and the prediction of even larger values, eH ^ 15toJ(^ lO^^G), at future LHC experiments 0, has also served to 
. inspire new studies on the effects of strong magnetic fields in QED and QCD. 

' On the other hand, the effects of strong magnetic fields in rclativistic theories is also important for planar condensed 
. ^ matter. The low-energy excitation quasiparticle spectra in systems like the pyrolitic graphites (HOPG) 0, ^-nd 
\l ■ graphene are characterized by a linear dispersion and then the dynamics of their charge carriers can be described 
I by a " relativistic" quantum field theory of massless fermions in (2-1-1) dimensions 0, [l^- Because the properties of 
. these systems can be modified in the presence of a strong magnetic field [llj , many works [l^, [I!] have been aimed 
at understanding the mechanisms behind the observed behaviors. 

Massless QED in the presence of a magnetic field exhibits a peculiar phenomenology. Due to the Landau quanti- 
zation of the fermion's transverse momentum in a magnetic field, the dynamics of the particles in the lowest Landau 
level (LLL) is (l-l-l)-dimensional. This dimensional reduction favors the formation of a chiral condensate, even at the 
d • weakest attractive coupling, because there is no energy gap between the infrared fermions in the LLL and the antipar- 
ticles in the Dirac sea. This phenomenon is known as the magnetic catalysis of chiral symmetry breaking (MCxSB). 
The MCxSB modifies the vacuum properties and induces dynamical parameters that depend on the applied field. 
This effect has been actively investigated for the last two decades In the original studies of the MCxSB 

[i31"[I3j the catalyzed chiral condensate was assumed to give rise only to a dynamical fermion mass. Recently, how- 
ever, it has become clear [ll] that besides the dynamically generated mass, the MCxSB inevitably produces also a 
dynamical anomalous magnetic moment (AMM), because this second parameter does not break any symmetry that 
has not already been broken by the chiral condensate and the magnetic field. The dynamical AMM leads, in turn, 
to a non-perturbative Lande g-factor and Bohr magneton proportional to the inverse of the dynamical mass. The 
induction of the AMM leads to a non-perturbative Zeeman effect [l^- An important aspect of the MCxSB is its 
universal character. It will occur in any relativistic theory of interactive massless fermions in a magnetic field. The 
MCxSB has been proposed as the mechanism explaining various effects in quasiplanar condensed matter systems . 

A drawback of the MCxSB phenomenon is that the dynamical parameters (mass and AMM) are extremely small 
even at relatively high fields. Hence, it would be important to have an independent way to experimentally detect this 
phenomenon. In a recent letter [l9| . we found that by measuring the induced electric polarization of the magnetized 
medium, compelling evidence in favor or against the existence of MCxSB can be obtained. 
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In this paper we find the non-perturbativc Eulcr-Hciscnberg Lagrangian of massless QED in a strong magnetic 
field, a theory in which the chiral symmetry is broken via the MCxSB mechanism and both the fermion mass and the 
AMM are dynamicahy generated [131 • it is known, the Euler-Heisenberg Lagrangian encodes the vacuum structure 
of QED in the presence of a constant electromagnetic field. This Lagrangian, first found by Euler and Heisenberg 
and then reformulated by Schwinger [2]| . is a low-energy effective Lagrangian describing the non- linear dynamics of 
the electromagnetic field. It is obtained from a single electron loop coupled to the external electromagnetic field. This 
formulation has been widely used to describe a variety of electromagnetic phenomena like photon splitting [22j , light 
scattering [2^, birefringence [l^], pair production [2^, etc. However, to the best of our knowledge, the low-energy 
Euler-Heisenberg Lagrangian has never been derived for a theory with dynamically generated parameters, which is 
the goal of this work. In the present case, the electron propagator in the loop will be the full propagator taken 
in the ladder approximation and in the presence of a strong magnetic field. Hence, our effective Euler-Heisenberg 
Lagrangian will be essentially nonperturbative because the electron self-energy that will be used to obtain it is found 
by a consistent resummation of an infinite number of rainbow diagrams. The nonperturbative origin of the effective 
action will be instrumental to show the paraelectric behavior of the strongly magnetized massless QED. 

The paper is organized as follows. In Sec. [Hi we present a brief review of the non-perturbative MCxSB phenomenon 
in massless QED. The dynamical parameters of the theory are found as solutions of the Schwingcr-Dyson equations 
(SDE) in the ladder approximation. These solutions will be later used to find the electric response of the magnetized 
medium to a weak electric field that is used as a probe. In Sec lIIIl we employ Ritus's eigenfunctions for constant 
and parallel electric and magnetic fields and the path integral formulation of Euclidean QED with a full fermion 
propagator, to obtain the non-perturbative Euler-Heisenberg action in a strong magnetic field. In Sec. IIVI we use 
the non-perturbative Euler-Heisenberg action to derive the electric susceptibility of the magnetized medium at strong 
magnetic field. Direct applications of the outcome of this paper to condense matter systems, like graphcne, are 
discussed in Scc.|Vl In Appendix [XI the Ep functions are obtained in Euclidean space for the case of parallel electric 
and magnetic fields. Finally, in Appendix B, we present an alternative approach, which starts from the effective (1-1-1)- 
dimensional theory of the electrons in the LLL, to obtain the Euler-Heisenberg Lagrangian of a strongly magnetized 
system. 



II. MxSB IN MASSLESS QED REVIEWED 



In this section, for the sake of completeness and understanding, we review the realization of the M^SB in massless 
QED, taking into account the dynamical generation of an AMM. In this approach, the electron Green's function 
depends non-perturbatively on the magnetic field, as well as on the dynamical mass, Af', and the induced AMM, T\ 
for each Landau level I. A more detailed description can be found in Ref. [isj . 

The Green's function that describes the motion of an electron in an arbitrary electromagnetic field satisfies in 
coordinate space the equation 

[j^U^^~j:{x,y)]G{x,y)=S\x-y) (1) 

where H^ ~ ^{df_i + ieA^) and 'E,{x,y) are the covariant derivative depending on the external field and the electron 
self-energy, respectively. 

In a constant magnetic field the fermion self-energy operator is a combination of the operators 7pH^, cr^yi*"^'^, 
(F^^^H^i)^ and (7^H'')^ [iSllBl- As all these operators commute with the later, the fermion self-energy operator will 
be diagonal in the basis spanned by the matrix-cigenfunctions Ep of (7pH^)^, 

(7 • n)2Ep = fEp . (2) 

For a magnetic field along the cca-direction, the eigenvalue p^, in the Landau-like gauge ~ (0, 0, iJxi, 0), is given 

by ~ {pq, 0, sgn{eH) y''2\eH\l, P3) with I = 0, 1,2, ... denoting the Landau Levels. From now on, we will consider 
in this section that sgn{eH) > 0. 

From the physical point of view, the Ep functions are the eigenfunctions of the asymptotic states of the charged 
particle in the presence of a constant magnetic field. The use of these eigenfunctions to diagonalize the self-energy 



operator in momentum space was originally done by Ritus [26j for the case of fermions. His method was later extended 
to charged spin-one fields in [23 . 

In the chiral representation, the Ep functions take the form 

Ep = ^ Ep,ix)Aia), (3) 

(T = ±l 
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where A(cr) = (/ + ia^'^^'^)/2 is the spin projector with a = ±1. The eigenfunctions Epa{x) are given by 



with Dn{p) the parabohc cyhnder functions of argument p = ^ 2\eH\{p2 / gH + x^) and index n = n{l, a) = I + ^-j^, 
and normahzation constant Nn = {ineHy^^ / \/n\. 



The Ep's are orthonormal [Ij 



d'^xEp{x)Ep, (x) = (27r)4,5(^) [p - p')Il{l) (5) 



and complete 



^-^Ep,ix)Ep,{y)^6Hx^y), (6) 

with notation E,p{x) = 7''Ej^7", and 

^ 5"'5{po-p'o)6{p,~p[)S{ps-p',), (7) 

n(/) = AisgnieH)) + A{-sgn{eH)){l - 5o0, (8) 

where n(Z) takes into account that for the lowest Landau level only one spin projection is allowed. Since we consider 
eH > 0, the spin projection at the LLL will be associated with A(+). 
The Ep functions satisfy the relations 

(7.n)Ep=Ep(7.p) (9) 

and 

+ Zip^)Ep(x) = Ep(Z||:^|| +Z^p^) (10) 



where p[[ ~ {po, 0, 0,^3) and pj^ ~ (0, 0, \/2eHl, 0), and, Z\\ and Z±_ are the longitudinal and transverse wave-function 
renormalization coefficients, respectively. 

Transforming to momentum space with the help of the Ep functions, it is easy to check, using Eqs. ((5|)- (fT0|) . that 
the fermion self-energy operator becomes diagonal in momentum space 

E(p,p') ^ |dWyEp(x)S(x,y)Ep,(y)-(2^)4r)(p-p')n(0S'(p) (11) 

with 

E'(p) = 4:^11 +zi:^^ + (M' + r')A(+) + (J\/' -r')A(-) (12) 

where M' and T' are the fermion dynamical mass and AMM for each LL, respectively. They have to be found self- 
consistently by solving the non-perturbative Schwinger-Dyson equations (SDE) of the theory. The presence of n(/) 
in Eq. ([TT|) ensures the separation of the LLL from the rest of the levels due to its lack of spin degeneracy. 
Applying the Ep transformation in ([T]), the full propagator in momentum space can be found as 

G\p,p') = j d^xd^yEp{x)G{x,y)Ep,{y) = {2nf^'-^\p - p')Ti{l)G\p) (13) 

where 

a.a=±l ' 

(14) 



4 




with 

AT 

N^{crT,aV\\) = aT^ - A/' - aV^\ 

D\aaT) = {M^f - {V^\ - aaT^f + {y[f 

V\\ = (i-4)iP||i 

y[ = {i-zi)\pA- (15) 

The non-perturbativc SDE for the self-energy can be solved using the ladder approximation, where it takes the 
form 

E(x, y) = ie]ff^G{x, yh^D^^ix - y). (16) 

Here, the full vertex is replaced by the free one, Df^,j{x — y) is the free photon propagator in the Feynman gauge, 
G{x,y) is the full fermion propagator (jl4p . and eji is the renormalized coupling constant at the scale \feH. In 
general, the truncation used for the resummation of the diagrams contributing to the SDEs in a given nonperturbative 
approximation can lead to gauge-dependent results. Therefore, care should be taken in choosing a gauge where the 
considered truncation is consistent and the physical parameters like the dynamical mass, the renormalized coupling, 
etc., are all gauge-independent. This is a well-known and old problem. For example, for the approximation used in 
the Abelian theory of fermions and vector mesons discussed in Rcf. [2^, the consistent gauge was the Landau gauge. 
In this special gauge all the vacuum diagrams with the Goldstone pole typical of theories with spontaneous symmetry 
breaking cancelled out, and as a consequence, the full vertex was replaced by the free one, so the truncation used 
there became reliable. In the presence of an external magnetic field, the consistent gauge in the ladder approximation 
is the Feynman gauge [l^ . One can show that in this gauge all the terms coming from loop corrections to the vertex 
cancel out and the full vertex reduces to the free one. The gauge-independence of the physical parameters in this case 
was verified in [H] by showing that in this special gauge the solution of the ladder SDE satisfies the Ward-Takahashi 
identities. Given that in the ladder approximation both the vertex and the photon propagator are replaced by the 
bare ones, the ultraviolet divergences that would appear in a pcrturbativc approach all cancel out, as long as the 
consistent gauge is used, so the photon field renormalization constant Z3, which also enters in the coupling constant 
rcnormalization, reduces to 1, leading to renormalized coupling and photon field at the characteristic scale of the 
infrared region where the SDE is solved, i.e. \/ eH (30| . Of course, any change in the approximation means taking 
a different truncation of the infinite series of diagrams and hence requires to find a new consistent gauge where the 
physical parameters result gauge-independent again. For example, going beyond the ladder approximation in the 
presence of a magnetic field, as in the so-called improved rainbow approximation that includes the polarization effects 
in the photon propagator, requires to use a non-covariant (Fcynman-like) gauge in which this approximation becomes 
reliable [U. 

One can show that the solution of is such that the mag netic field catalyzes the dynamical generation of a mass 
and an AMM even at the weakest attractive interaction [iSj ! This is the realization of the MCxSB [3] • Carrying 
out the Ep transformation in Eq. ([T5)) and taking into account Eqs. pT|) and (|13p we obtain 

lu.Y}i^m) - -IE ^ e-^i^e^("-""+-"-«')^J„„.(^.)^.".W.) 



(27r)4 g2 Vn!7i"!n'!n"! 
xA(a)7^A(a")n(r)G'"(F^)A(a")7,.A(a') 



(17) 



with n ~ n{l,a), n' = n{l',a'), n" — n{l",a"), n" = n{l",(j"), [a] meaning sum over ct, ct', a and a . In the above 
result we used 

D^%x -y)= I ^e'^^<^-y'>^ (18) 



(2714) q 



as well as \2£ 



dye-^'-%(y)7.Ep,(y) = (2^)^p)(y + - p)e"^e-^^^ ^ "^""'^^Vf . ' A(a)7.A(aO (19) 
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and 



J , ynml 



(20) 



where 



min{n,n') /i i— ln+n'-2m 

•■^ — ' miln — m)\(n — m)\ 

m—Q ^ ^ ^ ' 

and g_L = q^l2eH is the normahzed transverse momentum. 

Taking into account that in a strong magnetic field, the exponential factor in Eq. ([T7| serves as a cutoff for large 
transverse momentum, the main contribution to the fermion self-energy comes from the infrared region, g_L <C 1. 
It is then justified the approximation Jnn'{q±) — n\5nn'- Hence, the electron self-energy (jl7p in the strong-field 
approximation is given by 

5„,S'(p)n(0 = -z4EE / 7^^A(a)7''A(a")n(/")G'"(F^)A(a")7MA(<T')<S„„"<5w"„' ■ (22) 

I" [a] ^ ^ 

Now, using the relation = 6a-a-"Su" + 6-atT"Si"i+a- and performing the sum over I" and [a], we get 

E (p)n(O = -iej,j «^ <^ n(0 + 1 (23) 



Because the self-energy for a given Landau level I in Eq. p3p receives contributions from the fermion propagators 
depending on the Landau level I, together with the adjacent ones I — 1 and / -I- 1, the SDE's for all LLs form an infinite 
system of couple equations. However, the problem can be simplified taking into account that at strong magnetic 
fields, the leading contribution to each equation comes from the propagators with the lowest LLs, since the term IH 
acts as a suppression factor in the denominators of the fermion propagators with Z > 1. Under this approximation, 
one can find a consistent solution for each level [l^. Moreover, since all the solutions can be ultimately related to 
the LLL, which gives the main contribution, we have that all the dynamical quantities are determined by the LLL 
infrared dynamics. 

The above considerations imply that the SDE for the LLL reduces to 



S°(p)A(+) c ^^elA{+) j ^2 l'^\G\¥~q)ll (24) 

Substituting the LHS of ^ with Eq. we obtain that ^'fj = and 

where wc have used that 7ll^A^(p)7|j = 1 and introduced the notation for the rest energy £° = Al'^ + . Equation 
(j25p coincides with the one found in (3^ for the dynamical mass m^yn, because at the LLL, due to the lack of spin 
degeneracy, there is no way to separately determine the dynamical parameters and T", as they only enter in 
the LLL equation through the combination = -|- T". One can show [13 that the dynamical parameter 5°, is 

independent of the momentum in the infrared region <C f <C \eH\, and rapidly decreases, ~ as — > oo. 
This kind of behavior is characteristic of many dynamically induced solutions of SDEs. For example, it was also found 
for the dynamical mass in a theory with fcrmions coupled to a scalar field through a Yukawa vertex [l^ . 

The momentum dependence of the dynamical parameter S'^ indicates that the main contribution to the solution 
comes from the infrared region. This allows us to set p|| =0, cut the integral at the scale VeH, and assume that f° 
is momentum independent. After Wick rotating to Euclidean space, the LLL SDE becomes 

^ ' ^ (26) 



27r7 (27r)2 ^2 + 
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whose solution is given by 

£^ ~ V2i^e"v^ (27) 

where is the renormalized fine structure constant at the scale \feH. To keep our notation simple we did not write 
the subindex " R" in the electric charge " e" in the above expression, but from now on, in all the formulas related to the 
nonperturbative theory, " e" should be understood as the renormalized charge at the scale \feH . The non-perturbative 
character of the solution ((27)) is reflected in its dependence on a^. For higher LL's, each parameter, the dynamical 
mass and the AMM, can be found as functions of f ° [l^. They are all smaller than f °, reflecting the fact that the 
MCxSB is an infrared phenomenon essentially determined by the LLL contribution. We will use the solution (P7| to 
show that the chirally broken phase of massless QED in a magnetic field behaves as a paraelectric medium. 



III. NON-PERTURBATIVE EULER-HEISENBERG LAGRANGIAN FOR MASSLESS QED 



To investigate the electric response of massless QED in a strong magnetic field we are going to find the corresponding 
non-perturbative Euler-Heisenberg Lagrangian. The non-perturbative character of this effective Lagrangian comes 
from the fact that we are going to use the ladder full electron propagator considered in the previous section. Therefore, 
the mass and AMM parameters in our formulas should be understood as the dynamical solutions of the ladder SDEs. 
In all the derivations that follow, the ground state of the system is driven by the strong magnetic field, while the 
electric field just plays the role of a weak probe. 

The non-perturbative Euler-Heisenberg action <i>ey / \A\ can be found from the following path integral 



(28) 



where 



1 



(29) 



is the QED action in the presence of an external electromagnetic potential A^. Notice that it includes the electron 
self-energy T,{x), which in our case is the solution of Eq. (jl6p and thus a function of the dynamical parameters 
characterizing the phase with chiral symmetry breaking. Therefore, the effective action contains the correction to 
the classical QED action that comes from the resummation of the infinite number of rainbow diagrams entering in 
the ladder approximation for T,{x). To obtain an explicit expression for ^eff, we take (|29p as the initial action of 
an effective theory on which the gauge field fluctuations are neglected, and hence the electromagnetic flcld is merely 
reduced to a background field. 

In order to simplify the calculations, it is common to expand the effective action in powers of h. In the presence 
of an external electromagnetic field, the effective action can be expanded up to linear terms in h, where, since we are 
dealing with a non-perturbative theory, all the internal fermion lines are taken as the full fermion propagator. Then, 



eff 



$(0) +$(1) 



with 



d^x F^^F., 



(30) 



(31) 



and 



(32) 



The Fcynman diagram associated with (|32|) consists in principle of a fermion bubble with infinite insertions of photon 
fields [32|, which in this approach corresponds to external field lines. To investigate the linear electric response of the 
medium in a strong magnetic field, we are going to consider that in addition to the strong magnetic field, there is a 
weak but nonzero electric field. This last one can serve as a probe to explore the electric susceptibility of the system 
described by the action ([5^ . In Ref. [l^, we proved that the confinement of the electrons to the LLL by the strong 
magnetic field leads to an anisotropic electric susceptibility. In the plane transverse to the magnetic field direction 
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the susceptibility was just the same as in vacuum. Here we are interested in exploring the electric susceptibility in 
the direction along the magnetic field. With that aim we can consider constant parallel electric and magnetic fields 
in the cca-direction. In the Coulomb gauge they are described by the photon field = {—Ex3,0, Hxi,0). 

Performing a Wick rotation {xq — > —ix^, 70 — 5> — 474, and E — > iE) in Eq. ([5^ to go to Euclidean variables we 
obtain 



where for the sake of simplicity, we keep the same notation for the 7 matrices, etc, but it is understood from now on 
that they are all Wick-rotated. 

In order to transform (|33p to momentum space, it is convenient to use Ritus's approach and transform the spinors 
with the eigenfunctions Ep of the parallel fields case. 



ip{x) 



4^ i2ny 



Ep(x)^(p), ^(a;)=^^^b)Ep(x). (34) 

Details of the calculation of the Ep functions for constant and parallel electric and magnetic fields in Euclidean 
space can be found in Appendix [K\ There, we also verify some important properties satisfied by these functions. 
Next, using Eqs. (jA6p . (|A30p and (jA35|) . we can rewrite the Euclidean QED action in momentum space 

as 

/f-f\d'^vd^v' — — 
d^^E^-^^^i^{p)Ep{x)[^j^U^. - j:{x)]Ep,{x)i;{p') 

= %^iip)^''m"{i)[-i,-p, m]m (35) 

Here we took advantage of the fact that the self-energy S(p) is diagonal in the basis spanned by the Ep functions. 
We assume that the self-energy has the same structure as in the case with a pure magnetic field. This is a reasonable 
assumption, as we are interested in the situation of a strong magnetic field, but a very weak electric field, so the ground 
state of the system should not be affected by the presence of the electric field. Accordingly, the self-energy structure 



is still given by Eq. (fT2|) . but with the components of the momentum p defined by p|i = {sgn{eE)^ 2\eE\l, 0, 0, 0) and 

pj^ = iO,0,sgn{eH)^/2\ni\l,0). 

Carrying out the functional integrals in Eq. ([33|) . the Euclidean 1-loop effective action turns out to be of the form 

<^'e[A] = g ^ In [V'ip) + (M' + aT^] (l i^<5,„) (2 - 5 J (36) 



/,;=o 



with = (1 4- Z^^^){sgn{eE)^2\eE\l,0, {1 + Z^^_)sgn{eH)^2\eH\l,Q), (j)\,[A] = {Li)-'^<^\.[A], and L4 denoting the 
four-dimensional volume in Euclidean space. We call attention to the convenience of using Ritus's formalism to 
incorporate the non-pcrturbative effects coming from the self-energy E(p) in momentum space in a straightforward 
and simple way. 

As discussed in Sec. II, the ground state of the magnetically catalyzed system is mainly determined by the infrared 
dynamics of the fermions lying in the LLL. Out of this region, the LLL solution and with it, the rest of the 
dynamical parameters, quickly go to zero. Hence, it is consistent to limit the sum in / to the region where the chiral 
symmetry is dynamically broken, take the dynamical parameters as momentum-independent in that region, and just 
keep the leading contribution given by the LLL term (Z = 0, tr = +) in ((36|). Using the Schwinger's proper time 
representation for the logarithm and summing in 7, we find 

1 /-oo ^„ „ 2 

{<i,\)[A] = -e-^'^^ ^ \eH\\eE\coth{\eE\s) (37) 

StT^ Jl/\eH\ S 



where is given by Eq. (|27p . Notice that the consistency of the LLL approximation requires to introduce an ultraviolet 
cutoff l/|eiJ| in the s integration. This is in agreement with the fact, previously stressed, that the nonperturbative 
MCxSM is essentially an infrared phenomenon that takes place for energies below the scale \feH. For energies larger 
than this scale the dynamically induced rest energy goes to zero [30| . and the expression p7l) is not valid. 
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Let us make contact here with the usual Euler-Heiscnberg Lagrangian of pcrturbative massive QED in external 
electromagnetic fields. It can be obtained from p6p by replacing the dynamical parameters by the electron mass (i.e. 
A/' + (tT' m), introducing the Schwinger's proper time, and subtracting the vacuum contribution to regularize the 
divergence at zero fields, 

{^eYq^dIA] = [\eEs\\eHs\coth{\eE\s)cothi\eH\s)-l] (38) 

The effective action ([55)1 still has a field-dependent ultraviolet logarithmic divergence, which becomes evident in the 
expansion of the hyperbolic cotangents in the squared parenthesis. Following Ref. [2]| . the renormalization procedure 
is achieved by renormalizing the charge and fields and subtracting in the effective action the counter term, 

^,yrr^A^ . ^^"^ ,39, 

where the charge and fields appearing in Eqs. ((38)) - (p9)) should be understood as the renormalized ones. 

Subtracting p9p from psp . and doing the analytical continuation E — > —iE, we recover the well-known Euler- 
Heisenberg Lagrangian of massive QED in Minkowski space [s^ . 

The path integral formulation of the effective action for massive and masslcss QED in the presence of a constant 
electromagnetic field has been previously found in Refs. [H, i35i] in terms of the spectra of the Euclidean Dirac's 
operator and the generalized ^-function, but within a perturbative approach. A pcrturbative approach docs not allow 
to incorporate the phenomenon of MCxSB. However, the MCxSB is unavoidable in the massless case, because once a 
magnetic field is present, no critical magnetic strength is required for the dynamical breaking of the chiral symmetry 
to occur and for the mass and AMM to be generated. 



IV. PARAELECTRICITY 



With the help of Eq. (|37p we can now show the paramagnetic behavior of massless QED under an applied magnetic 
field. First, notice that the electric polarization of the system can be obtained as minus the derivative of the effective 
action, or electromagnetic free-energy density $ with respect to an applied electric field. For a weak electric field E, 
the free-energy density can be expanded in powers of E as 



(40) 



In a strongly magnetized medium, the coefficients $07 Xi etc., may in principle depend on the magnetic field. The 
susceptibility coefficient r] is different from zero for ferroelectric materials [36j . In magnetized QED it is zero, because 
the second term in the RHS of Eq. violates parity, a symmetry that is not broken neither in massive QED nor 
in the chirally broken phase of massless QED. Then, the coefficient x characterizes the lowest order of the system 
dielectric response. It accounts for the electric polarization P = xeE = ^ induced by an externally applied electric 
field. This term does not break any additional symmetry of the free energy 

Assuming that E ^ (^ °)^, we can perform a Taylor expansion of (jeii'ls) cothdei^ls) in power of |ei?|/(£°)^ in ()37p . 
Keeping just the leading term in the expansion, the effective action can be approximated by 



\eH\ 



ds 



-s{£° 



\eE\s) 



21 



\eH\{\eE\f f'^'" 



247r2 



ds 



l/eH 



oiR \eH\ 



E' 



(41) 



Note that the Feynman diagram associated with this equation corresponds to the photon polarization operator with 
internal fermion lines of full fermion propagators, two external legs of "electric" photons, and an infinite number of 
insertions of "magnetic" photons. A st raig htforward calculation of the polarization-operator in the static limit 
corroborates the above observation [l^ l37l |38| . 

To regain the effective action in Minkowski space we should make the analytic continuation, E — >■ —iE, in (|41l) . 
Comparing the obtained action with Eq. pO)) . we immediately find that for the magnetically catalyzed QED the electric 
susceptibility in the direction of the magnetic field is 



Xe 



OCR |eg| 

3^ (£0)2 



(42) 
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Using ([27]) in ([42]) we obtain 



Xe = T^eV "R , (43) 

OTT 



in agreement with the rcsuh found in the infrared hmit from the photon polarization operator in Ref. [l9[ . The non- 
pcrturbative dependence of the susceptibility on the renormalized fine-structure constant an accounts for the large 
electric response along the field direction of the medium with MCxSB. The susceptibility in the transverse direction 
was found in [l9| . It resulted to be zero as in vacuum. Therefore, at strong magnetic field the system with MCxSB 
displays a noticeable anisotropy in its electric response. 

It is worth to point out here the difference between the electric responses of the MCxSB system and of the 
regular massive QED in a strong magnetic field. In the massive QED case, the Euler-Heisenberg action is found 
perturbatively. The one-loop diagrams contributing to the perturbative action contain a logarithmic divergence that 
needs to be regularized, as pointed out in the previous section. This ultraviolet divergence affects the charge and 
field renormalization and leads to a logarithmic contribution ^ ~ H^)\xi{eH /m?) to the action, which in turn 

affects the susceptibility in the parallel and transverse directions [3^. In this situation, the strong field limit can be 
taken only after the divergences of the theory are eliminated. In the nonperturbative case studied in the present work, 
the ultraviolet logarithmic divergence is absent because the dynamics of the fermions that leads to the magnetically 
catalyzed, chirally broken phase is essentially infrared {p^ « eH), so the nonperturbative Euler-Heisenberg action 
has a dynamical cutoff at the scale VeH that limits it to the infrared region where the dynamical parameters are 
momentum-independent and nonzero. Besides, in the ladder approximation used to obtain the physical parameters, 
the full photon propagator and vertex reduce to the free ones, so any potentially dangerous infrared divergence is also 
cancelled out and the electromagnetic field and coupling constant are automatically renormalized at the scale \feH . 
Hence, in the nonperturbative case the transverse susceptibility remains zero as in the vacuum. 

It should come as no surprise the different electric behaviors in regular massive QED and in massless QED with 
MCxSB. In massive QED, the ground state is driven by unpaired electrons, while in the nonperturbative case, the 
electrons pair with positrons forming tiny electric dipoles. As the electrons forming the dipoles are all in the LLL, 
their dynamics is (l-f-l)- dimensional (their dispersion only depends on the parallel momentum) due to the well-known 
dimensional reduction that occurs at the LLL. Hence, the electric dipoles can be polarized by any small electric field 
to produce a strong polarization only in the parallel direction. The role of the magnetic field here is to induce the 
dipole moments, while the role of the electric field is to polarize them. 

The lack of explicit magnetic field dependence in the susceptibility (j43p is a unique feature of the paraelectricity 
in massless QED with MCxSB. This propert y is not found in other strongly magnetized systems like massive QED, 
nuclear matter, color superconductivity, etc. [40[. Even in the case of MCxSB in the context of QCD in the presence 
of a strong magnetic field, the c/iro?7io-susceptibility remains an explicit function of the magnetic field through the 



running of the strong coupling [41 



V. CONCLUSION 



In this paper we obtained the non-perturbative Euler-Heisenberg Lagrangian for massless QED with MCxSB. 
The magnetic field catalyzes the breaking of chiral symmetry through the generation of the dynamical parameters 
that are solutions of the non-perturbative Schwinger-Dyson equations. Since the dynamics of the electron-positron 
pairing leading to the breaking of the chiral symmetry is mainly induced in the region of momenta much smaller 
than V eiJ, the magnetic field introduces a dynamical ultraviolet cutoff in the theory that also enters in the non- 
perturbative Euler-Heisenberg action and provides the scale at which the coupling is renormalized in the calculations. 
The electric response of the medium with MCxSB, characterized by the electric susceptibility, can be obtained from 
the Euler-Heisenberg Lagrangian. The non-perturbative character of the calculation is reflected in the dependence 
of the susceptibility on the coupling constant. A remarkably feature of the electric response is its anisotropy. In 
the direction parallel to the magnetic field, the susceptibility is large and independent of the magnetic-field strength, 
while in the direction transverse to the field it is zero. Given that the amplitude of the dynamical parameters are 
typically too small to be observable in the experiment, the large electric susceptibility is the best candidate available 
to test the realization of the MCxSB phenomenon in a given physical system. 

An important implication of this result is that the chirally broken phase exhibits strong paraelectricity, a property 
found in certain condensed matter systems like quantum paraelectric (QP) materials [l^] and transition-metal-oxides 
(TMO) [i^. In those materials, unaligned electric dipoles are aligned in an external electric field, producing a high 
electric susceptibility, often exceeding 10^. In QP materials the large electric susceptibility is temperature-independent 
below certain critical temperature, a property attributed to a quantum phase transition [4^ . An interesting question 
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to explore in the future is whether the strong susceptibihty found here within a (3+l)-diniensional theory is also 
present in quasiplanar condensed matter systems as bilayer graphene. It is known, that the band structure of bilayer 
graphene can be controlled by an applied electric field perpendicular to the layers' plane. The electric field creates an 
electronic gap between the valence and conduction bands with energy values that varies from zero to mid- infrared (4^ , 
depending on the field strength. Under a very weak electric field the gap is practically zero and the spectrum is Dirac- 
like. Even though this is a very peculiar (3+l)-D system, only formed by two layers, one could attempt to model it 
with a (3-|-l)-D theory of interactive masslcss fcrmions. Due to the universality of the MCxSB, we expect that the 
application of a strong magnetic field parallel to the weak electric one will trigger the generation of a dynamical energy 
gap. Under these conditions, one would expect that detecting a very large electric susceptibility in the direction of 
the applied fields would signal the realization of the MCxSB phenomenon. 

Another interesting direction worth to be explored in the future within the framework of the non-perturbative 
Euler-Heisenberg Lagrangian is in connection with higher order non-linear effects. Recently [45| . the study of higher 
order nonlinear effects in the context of the perturbative QED Euler-Heisenberg Lagrangian indicated that a purely 
magnetic moment placed in an external quasistatic electric field can lead to the induction of an electric moment and 
viceversa. It would be interesting to investigate if a similar nonlinear effect occurs in the case of the magnetically 
catalyzed system. 
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Appendix A: Ritus Ep functions for parallel electric and magnetic fields in Euclidean space 



In this Appendix, we find the matrix eigenfunctions Ep for charged fermion in the presence of constant electric and 
magnetic fields. They can be used to diagonalize the full fermion propagator in momentum space. This approach has 
been extensively used in the literature for the case of a constant magnetic field, where the Ep functions play the role, 
in a ma gnet ic field, of the plane waves (Fourier transform) at zero field. This is the essence of the so-called Ritus' 
method [2g. Here, we will consider this approach for the case of constant and parallel electric and magnetic fields, 

X i? = in Euclidean space, where these fields enter in a very symmetric way. 

In the presence of electric and magnetic fields, the self-energy operator is a combination of four scalar operators: 
7^n^, c7,,,F^^ (F^'^H^)2, and jsF'^'Ff,,, with F^, = d^A, ~ d,A,„ F^"^ = ie^,„,F"^ a^, = and 
H^ = {idf^ — eAf^). Since all these scalars commute with (7^H'')^, the self-energy, and hence the fermion propagator, 
will be diagonal in the basis spanned by the eigenfunctions of (7^H'')^. 

Without loss of generality, we assume that the fields point along the positive xs-axis. We fix the gauge, = 

(— Fzs, 0, iJxi, 0), for the external field, and work in the metric g^i^ = (1, — 1 ). Taking into account that [np,n^] = 
—ieFfj^^, one can write the operator (7^H'')^ as 



t"3, 

= {ido + eExsf - {id^f - {id2 - eHxif - {id^f - e (^£75X3 - Sgi?) (Al) 
Working in the chiral representation of the gamma matrices 

^"=(-1^0^)' ^'=(-^0)' ^■5 = (o-i) (^2) 

on which both S3 = 17172 and 75 = 170717273 are diagonal, the eigenfunctions of the square Dirac operator eigenvalue 
problem 

(7^H^)V=pV (A3) 
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can be written as 

= Ep^„^r{x)Vcr,r, (A4) 

where the functional part Ep^„^r{x) of the eigenfunction has to be found from the eigenvalue equation (jA3p . while 
the spinors v^r are chosen as the eigenvectors of the spin S3 and electric dipole 75S3 operators, i.e., 'E^Va-.r = o't'<T,r, 
"f^HsVa^r — rvcr^r, with (7 = ±1 and r = ±1. 

Since we will use the Ep functions in the framework of the path-integral formulation, we just need to find them in 
Euclidean space. With this aim, we perform a Wick rotation xq — > —ix^, so that the Euclidean eigenvalue equation 
becomes 

[(-^^4 + eExsf + {-tdif + {-td2 - eHxif + {-id^f -e{rE + crH)] Ep,^,r{x) = p^Ep,^,r{x) (A5) 
To solve Eq. (|A5p we separate in longitudinal and perpendicular variables via 

Ep,,,r{x) = E^^,{x\\)E^^,{x^) (A6) 

where 

£;^^,(x|i)EEe'^'-^X^.(^3) (A7) 



E;^,,{x^) EE e^-P^^-xlM) (A8) 
and substitute Eq. (jA6p into Eq. (|A5|) to obtain an eigenvalue equation for each x function. The one for ^[x\) is 

{-dl + b2 - eHx^f - (pi + eHo)\xlAxi) = (A9) 

where the eigenvalues 'p\ = ^ ' ^'^'^ P\\ ^'"'^ constants usually introduced in the separation of variable method 
that has to be determined from the eigenvalue equations. Eq. (jA9p is the same as the one found in the case of a 
system with just a magnetic field (29j . Its solutions are the parabolic cylinder functions Dn„{pH) with argument 
Ph = ^2\eH\ — xx) and index 

UH = nH{(T,l) = l + sgn{eH)^-]^, 71^ = 0,1,2,... (AlO) 

The non-negative integer I is the Landau level, which as known characterizes the quantization of the transverse 
momentum in a magnetic field. Notice that in the LLL, I ~ 0, only one spin projection is allowed, that is, a = +1 if 
sgn{eH) > 0, or cr = — 1 if sgn{eH) < 0. The eigenvalue satisfies = — = 2\eH\l. 

The equation for Xp,r{xi) is 

-^3 + iP4 + eExsf - (pfi + eEr) 

which, as Eq. (|A9|) . is also the parabolic cylinder equation and hence has solutions of the form Xpri^s) = DnEips) 
with argument pE = \/2\eE\ {x-^ + and index 

~ ~ r 1 

HE = nEir,l) = l + sgn{eE)-~ = 0, 1, 2, ... (A12) 

On the other hand, the non-negative integer I enters in a similar way to the Landau level, but it corresponds to 
the quantization of the longitudinal momentum in the presence of an electric field in Euclidean space. Therefore, in 
Euclidean space, there is a symmetry between the electric and magnetic field sectors of the solution. Notice that r 
represents the projection in the electric field direction of an "intrinsic" electric dipole moment. Only one projection 
of r is allowed at Z = 0. _ 

From Eq. (jAlip one finds p| = 2|ei?|Z, so solving for we obtain for the eigenvalue in (jASp . p^ = Pj_ + p| = 

2|e_ff |Z + 2\eE\l. The corresponding normalized eigenfunctions are 

Ep,aA^) = Nr,„^P^'''Dn„{pH)N^Ee'P''''D^,{pE) (A13) 



(All) 
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with normalization constant Nn^ = (47r|e-ff |) * / V'^h'i a-nd Nn^ found by replacing by i? in Nn^^. 
Introducing the spin and electric dipole projectors 



A^(a) = i(l + aS]3), A^(r)==i(l + r75S3) (AM) 



which satisfy 



A^(+)A^(-) = 0, A^(+)A^(^) = (A15) 

A^(+) + A^(-) - 1, A^(+) + A^(-) = 1 (A16) 

[A^(a),A^(r)] =0 (A17) 
The Euclidean Ep functions of the problem with parallel constant electric and magnetic fields can be defined as 

E, =E^(x||)E^(.Tx) (A18) 

with 

E^(x||) ^ E^.r{H)A^ir) - e^^^^^ E X^.(2;3)A^(r) (A19) 

r=±l r=±l 

Ef (x^) ^ J2 E^A^i_)A"ia) = e^P^-^ ^ X^Ax,)A" (a) (A20) 

(T=±l cr=±l 

It is known that spinors in Euclidean space obey 

{ipix),ipiy)} = mx)My)} = Mx)My)} = o (A21) 

where the last relation implies that ip is not necessarily obtained as the product (^^74. Then, (p and result in two 
totally independent functions that should be found from two corresponding independent Euclidean equations. This 
is one of the main novelty of Euclidean fermion field theory [i^ . 

Hence, to obtain the Euclidean Ep function, since it is associated with the spinor Tp = 'Pa,rEp^cr,r, we need to find 
Ep^a,r as the solution of the Euclidean eigenvalue equation 

[(^94 + eEx3)^ + {idif + {182 - eHxif + {id^f ^e{rE + aH)] Ep^,^r{x) = p^Ep^^A^), (A22) 

Once wc performed the procedure shown in Eqs. (jA5[) - (jA20p . we arrive at 

Ep = Ep{x\\)E"{x±) (A23) 

with 

E;(x||) ^ J2 ^^.(^ll)A^(r) = e-'^^-^ ^ X^.(^3)A^(r) (A24) 

r=±l r=±l 

e"{x^) ^ J2 EpA^^)A''{a) = e-^^-^^ ^ X^.(xi)A^(a) (A25) 

<J = ±1 <T=±1 

To finish, it is worth to mention that the solutions of (jA3|) in Minkowski space has been worked out in [2^, 113, . 
They are widely used in the scattering matrix method. In Minkowski space, the charged particles in the presence 
of a uniform electric field do not have bound states. However, after the Wick rotation to Euclidean variables, the 
dispersion by a barrier problem is transformed into that of a particle in a potential well. The solutions to the field 
equation in this case are instanton-type, that is, classical solutions in Euclidean time [iol. [soj. 
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2. Properties of the Euclidean Ep functions 

Starting from (jA18[) and (|A23[) . it is easy to check the orthogonaUty and completeness relations of the Ep functions. 
One can explore the orthogonality by separating the integration in parallel and transverse coordinates 

d^xEp{x)E.p'{x) = J d^xii Ep (a;||)E^(a;||) J d'^x± {x±)E^, {x±) (A26) 

Let us focus our attention on one of the integrals in the RHS of Eq. (|A26p . Let's say the one associated with the 
electric field. 

yd^xii Ep(a:||)E^(x||) = J2 N^ENnE'^''{r)A''ir') J dx^dx3e'^P^'-P-^^-D^,{pE)D„,^{p'E) (A27) 

r,r' — ±1 

where ps = \/'2\eE\{x3 +p4/eE), p'^ = ■\/2|ei?|(x3 + P4' /eE). The integration in X4 produces a delta function 
d{pi — Pi), which implies that pE = p'^- 

Changing variables from X3 to pE and using the orthogonality of the parabolic cylinder functions 

j dpEDj,j^{pE)Dn'^{pE) ^ V27rnE^-SnEn'j^ (A28) 

we obtain 



d^xii Ep(x||)E^(x||) - {2n)S{p4' - Pi) ^ Nn.N^.-Srr' ^— 



nElSn^n'^^ir) 



= {2TT)^Sii,S{pi' -pi)n''{l) (A29) 

where n^(Z) — {sgn{eE)) + {—sgn{eE)){l — (5;q) takes into account the condition sgn{eE)r = 1 if / = 0. Note 
that a similar equation is obtained in the external magnetic field case. This similarity between the magnetic and the 
electric fields, is a consequence of the duality between the electric and the magnetic fields which is evident from the 
equations of motion in Euclidean space. 

A similar procedure can be carried out for the magnetic part in Eq. (|A26[) . Thus, the Ep functions satisfy the 
orthogonality condition given by 

d^x Ep{x)Ep,{x) = {2n)*5u,Su'5ipi' ~pi)Sip'2-p2)Il''il)Il"il) (A30) 

where the projectors n^(Z) and n^(/) guarantee that in the quantum states with I = and/or I = 0, respectively, 
the particle can have only one spin and/or dipole moment projection. 

The completeness relation can be easily proved by using the orthogonality condition of the Ep functions, as follows. 
Let us multiply the RHS of Eq. (jA30p by Ep'(y) and then integrate over p', 

J d^xEpix)^-0^Ep,ix)Ep,iy) - ^dV%<5„,<5(p4' - P4)<5(p^ - P2)n^(0n^(0Ep' (y) 

= n^([)n«(Olp(y) = Ep(y). (A31) 
Thus, from the above equation, we can see that the Ep functions satisfy 

'^\Ep,ix)Ep,{y)^d\x~y), (A32) 



'(27r) 

which is the completeness relation of the Ep functions. 

From the orthogonality condition of the Ep functions we have that the condensate does not depend on the repre- 



sentation. 



d^x i;ix)ijix) = ^^iip)^{p) (A33) 

with l^gr - Y.rj=o I-oc dPidP2- 

Relation (|A33P follows automatically from the Ep-transformation of the wave functions of charged particles 

^(^) = ^ (04Ep(^)V'(rt, = ^^mM^) (A34) 

together with the orthogonality condition (|A30p . 
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3. Generalized Momentum for the Parallel-Field Configuration 

To find the generalize momentum of the charged particles under a parallel constant electromagnetic field we 
should solve the equation 

yiEp = Ep7 • p (A35) 

Taking into account Eqs. (|A6|) - (|A20p . we have that the particle dynamics along transverse and longitudinal coor- 
dinates are decoupled. Then, the LHS of Eq. (|A35|) can be written as 

mp = (P" + n^)^p^p - (p"E^)Ef +E^(p^Ef ) (A36) 

where we took into account ()A6[) . and that Ilj] = {—id^ + eEx^, 0, 0, —id^), ~ (0, — —182 ~ eHxi,0). 
Now we can separately solve the equations 

pliE^ = E/(7ll.|5||) (A37) 
P^E« = Ef(7^.p^) (A38) 

Let us focus our attention on the LHS of Eq. (|A37p . Without lost of generality, let us assume that sgn{eE) > 0. 
Thus, using that 73 = — ^237574, Eq. (|A37|) can be written as 

pIlE^ = e^P^^^ { [aA^H + at A^(+)] 7^} ^ .(^3)A^(r) (A39) 

r=±l 

where a = {p^ + eEx^) + and = [p^ + eEx^) — are lower and raiser ladder operators, respectively. 
Taking into account that A^(r)74 = 74A^(— r) together with 



atx^_ = V2|ei?|fx^+ (A40) 



ax^^+ = ^J2\eE\l x^- , (A41) 



which are easily obtained by using aDj{7]) — ly/2\eE\Di_-^{ri) and a'' DiZiiv) — A/2|e£'|D;-(r/), we can rewrite Eq. 
(|X39)) as 



with 



A^i~)^2\eE\l xt^ + A^{+)^J2\eE\l xl 



7^=E^(7ll-P||) (A42) 



pll = (V2[|eii;|, 0,0,0). (A43) 

where we used the fact that for / > 1. n^(^) commutes with 74, and, when I = the RHS of Eq. (|A42p is identically 
zero. 

A similar procedure can be followed with Eq. (jA38p to obtain 

p^ = (0,0,V2Z|ei?|,0). (A44) 
Substituting these results back in (|A36|) we have 

pEp = Ef Ef (7!! • + E^E^(7^ • p^) = Ep7 • p (A45) 
where the generalized momentum, in the general case, is given by 



Pp=P^l+Pp = isgn{eE)^2l\eE\,0,sgnieH)^2l\eH\,0). (A46) 
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Appendix B: Non-perturbative Euler-Heisenberg Lagrangian for electrons in (l+l)-dimensions in a weak 

electric field. 



As we already mentioned, in the presence of a strong magnetic field charged fermions suffer a dimensional reduction 
from (3+1) to (1+1) dimensions. This is because all fermions are confined to the LLL. In this section we shall follow 
an alternative approach to find the Euler-Heisenberg Lagrangian of the studied system starting from an effective 
theory of fermions in the reduced (1+1) dimensions under a weak electric field. 

Let us start from the Euclidean action (I55t 



Because the factor (l) emphasizes that fermions in the LLL have only two degrees of freedom compared to the 
four degrees of freedom of those in higher LL's, let's separate the contribution of the LLL from the rest 

s = So + ]^ ^ (04^(f )n''(0h7MPM - mWp) (B2) 

where 

= E / ^#^(p)n^(z)A(+)[-7^i5^ - m]m (B3) 

is the action of fermions in the LLL with Ie = l/VeE and Ih = 1/VeH being the magnetic and electric characteristic 
length scales, respectively. Note that n^(^) plays the same role as its magnetic counterpart n^(^) separating the 
dynamics of fermions in the Lowest Electric Level from the rest. However, because we are interested in a scenario in 
which the electric field is weak, we keep all electric levels in So. 

Using the projectors of Eq. (|A14|) . the spinor field can be decomposed as follows 

= A^(+)A^(+)^> , Vi+)=A^(-)A^(+)^, Vi") = A^(-)A^(-)^> , = A^(+)A^(-)V , (B4) 

V^^+^=^A^(-)A«(+) , ^i+^=M^(+)A«(+) , ^^-^=^A^(+)A^(-) , ^i-^-^A^(-)A«(-) . (B5) 

where the supraindices (±) denote the spin-up (+) and spin-down (-) projections, while the subindices (R/L) are 
labeling the right (R) and left (L) chirality projections. 

In lEl) and dH]) it was used that A^(cr)A^(r) = A"{(T)P{ra) with 

P{ar) = \{l + r^775) , (B6) 

the usual chiral projector: ~ P{+)'4' and ipL = P{'—)4'- 
It can be easily checked that 

^ = + V-i^^ + ^k"^ + i'V ■ (B7) 
Then, using Eq. (|B4p and Eq. (|B5p . we rewrite the effective action for the LLL of (|B3p as 



5*0 — y 

[=0 



- f " [?^k+^(p)H^(04+)(p) + V^i+^(p)H^(0^4+)(p)] } . (B8) 



where S(p) = A(+)£0. 

Writing the four Dirac spinor in the LLL as ifF = {ipi,ip2, ip3,ip4), Eq- (|B8[I reduces to 



^Elnf^J (2^)4 \-iJ2\eE\l -£° J 



where n(Z) ~ i(l + <735j^) with (73 the Pauli's matrix. 
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Introducing the Dirac gamma matrices in the (1+1) Euclidean space 

74 = i<Ji , 71 = -10-2 and 75 = 0-3 = 17174 , (BIO) 

which satisfy {7,j,7;y} — —2(5^^ with 5^i, — diag{\, 1), and the 1+1 spinor field in the LLL as 

then we rewrite Eq. (jB9p as 

^" = E / ^^LLLn^(0 {~1,P, - £') i'LLL (B12) 

/=0 



where = {^2\eE\l,Q). 

Note that the above equation in the limit — ;> describes a theory of fermions in the LLL of the MCxSB phase 
[l8| . This statement becomes evident once we make the replacement {\/Ie) Sr~^ / '^Pi^ ^^id — > p|i = (^4,753). 
Taking into account Eq. (|B12[) . the Euclidean 1-loop effective action <^\, for electrons in the LLL has the form 



e*- = ][D{i^LLLm{i^LLL}]e'^ 

= Vet [n^(0(-7^p^ - f ")] . (B13) 
We can rewrite the effective action in Eq. ()B13P as 

1=0 

where (f)\^ = (Z_e///)^^$]j and we took into account the definition of n^(/). Once we integrate out and p2^ we get 

\eH\ \eE\ . . f°° ds 



\eH\ ds 



e 



f rl 

/ ^e-''^^°^\\eE\s)cot\v{eEs) (B15) 



2(2^)24 s2 

where we used the ultraviolet cutoff Pj^ in s that accounts for the LLL dominance of the pairing dynamics in this 
theory. Note that Eq. (jB15|) coincides with Eq. ((37)) as expected. 
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